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Abstract 

We address the question of understanding the production of J /iff particles regarding the global underlying event in 
proton-proton collisions. To do so, we propose to look at a new observable: the J/i// production as a function of the 
charged particles multiplicity of the event. We demonstrate the interest for an experimental measurement by varying 
the model of multiple interactions in the PYTHIA generator. 
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1. Motivation 

Understanding the production of the J /if/ particle in 
proton-proton collisions is a complex issue ([llSl- It 
becomes even more complex when one considers the 
production of a J/i// in a complete event containing in- 
terplaying soft and hard components produced simulta- 
neously. In this first section, we provide motivations for 
studying quarkonia production regarding the complete 
proton+proton event. We emphasize that proton-proton 
collisions are not trivial, especially in the regime of high 
multiplicity events. 

Very few observables have been studied as a func- 
tion of the charged particles multiplicity. The CDF col- 
laboration conducted such a study for the mean trans- 
verse momentum (pr), as illustrated on Fig [1] An im- 
portant raise of (pr) is seen as a function of multiplic- 
ity. This can result from several elementary multiple 
interactions occurring in parallel in growing multiplic- 
ity events. In these events, one can also expects non 
trivial phenomenon, such as collective effects, as men- 
tioned later. Thus, there is no theoretical clear picture 
for the explanation of {pj) versus multiplicity, but mul- 
tiple interactions clearly play an important role. An- 
other feature that helps to understand the importance 
of multiple interactions, is the fact that, for proton- 
proton collisions, the jet cross section computed in the 
framework of perturbative Quantum Chromodynamics 
(pQCD), becomes larger than the total cross section: 



there is thus more than one elementary collision in high 
multiplicity events ^ . 




Pythia hadron level : 
■TuneAnoMPI TuneAp^=0 
- ATLAS tune 



TuneA p^=1 .5 
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Figure 1: {pj) as a function of multiplicity. The data points (black 
squares) are compared to different event generator simulations from 
PYTHIA. In PYTHIA: improvement of the multiple interactions 
scheme leads to a better description of the CDF data (J. 

On Fig. [1] authors compared this observable with 
the event generator PYTHIA fsl). These event gener- 
ator is based on pQCD and assume the factorization 
scheme. At the time the CDF data. Run I, was re- 
leased, none of the models were able to describe the 
data i6ll- After improving its treatment of multiple in- 
teractions, PYTHIA finally describes better these data 
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with the tune A {pj_ - 1 .5) 101 ■ A good understanding of 
high multiplicity phenomena is thus important to com- 
pletely describe events with hard and soft components, 
and possible interactions between the two. It is to be 
noted that collective effects in proton-proton collisions 
can also contribute to explaining this observable, as for 
instance accounted in the EPOS event generator 101 ■ 

The only other multiplicity dependence study that we 
are aware of is from one of the first UAl paper where 
the authors studied the transverse energy Et as a func- 
tion of multiplicity Besides these CDF and UAl 
observations, and a recent CMS result that will be men- 
tioned hereafter, high multiplicity phenomena in proton- 
proton collisions are not studied. On the opposite, in 
heavy-ions collisions, studying observables as a func- 
tion of the (very high) multiplicity is a key concept and 
a very common tool to study the Quark Gluon Plasma 
(QGP). 

2. J/i/f production in complete events 

In standard proton-proton collisions, J /iff are pro- 
duced through hard processes. Hard processes can be 
understood in the framework of pQCD, where in the 
limit of high momentum transfer Q^, a hard compo- 
nent (2 — » 2 or 2 — > 1 processes) can be factorized 
out from the soft component, as predicted by the factor- 
ization theorem Jgl [l3l ■ J/ 4' ^re very particular among 
the hard probes, because they first imply the production 
of a massive charm quark- antiquark pair, followed by 
its binding into a charmonia (cc) state. Nowadays, the 
exact charmonia production scheme is not clear and is 
debated within the community (for reviews on this as- 
pect, see in |2l)- Furthermore, no models describe the 
underlying event associated to a J I if/ production, and 
one needs to be aware of this lack of description, when 
dealing with event generators. This problem was for in- 
stance raised in February 2010 during a workshop on 
quarkonia production at CERN iflTll . New observables 
are then needed to start understanding and constraining 
the interplay of soft and hard aspects of proton-proton 
collisions, the importance of multiple interactions in the 
production of hard processes and the possibility of new 
phenomena such as collective eff'ects. 

The simple new observable we propose to look at is 
the J I if/ yield versus the underlying event multiplicity. 
The basic idea is to compare the multiplicity of events 
containing a J/ip to the multiplicity of Minimum Bias 
events. The physics that we want to address here is the 
interplay between the hard {J /if/) and the soft (underly- 
ing) components in a complete event, the influence of 



one over the other, in the case of quarkonia production 
such as the J /if/. The question is: do the underlying 
event and the multiple interactions scheme matter in the 
production of hard processes, such as a J/if/1 Though 
we limit ourselves to the case of J /if/, we think the com- 
plete understanding of global events will only arise by 
studying as many quantities as possible, including vari- 
ous hard probes. 

3. First simple PYTHIA study 

To test the relevance of our new observable, we con- 
ducted a simple study at the event generator level. To 
do so, we needed an event generator that is able to pro- 
duce simultaneously the J /if/ and the corresponding un- 
derlying event. PYTHIA is one of the most developed 
event generator having the two functionalities. Other 
event generators have more sophisticated description of 
the soft component (HERWIG [H, EPOS fTj), but 
none of them include the production of quarkonia. An- 
other possibility could be the CASCADE event genera- 
tor lll4ll . this could be explored. 

We generate two uncorrected sets of PYTHIA 
events. In the first set, it is required to produce a J/if/ in 
each event (ISUB=86), through the color singlet model 
used by default in PYTHIA ||5|] . The second set is made 
of Minimum Bias (MB) events (MSEL=2). The mul- 
tiplicity distribution is determined for the two sets and 
their ratio {J/if/-^ MB) is computed, then normalized for 
the J/if/ and MB cross sections, estimated by PYTHIA. 

It is to be noted here, that our J /if/ are directly pro- 
duced and do not come from B,Xc or if/' decays. 

Fig. [2] shows this first result with a milhon gener- 
ated events for each set, for proton-proton collisions 
at 10 TeV. Here, Multiplicity stands for the number of 
charged particles with pr > 0.9 GeV/c and pseudora- 
pidity lower than 2.4, corresponding to the acceptance 
of the CMS tracker. In all the PYTHIA tunes we consid- 
ered, the J /if/ production mechanism is the default one, 
meaning the Color Singlet model, but one should also 
looked at different implementations of the production 
mechanism for other birth processes (such as the Color 
Octet mechanism). 

The four curves refer to diff'erent PYTHIA settings, 
mostly implying variation of the multiple interactions 
model. 

• The green closed circle symbols account for basic 
PYTHIA, version 6.2, with all parameters set as 
default (MSTP(82)=1). In this simple model aU 
interactions are supposed to be independent. The 
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number of parton-parton interactions in a proton- 
proton collision follows a Poissonian law, with the 
mean number of interactions h - cjiard/iTineiastic- 
Here, the hard process is the first one, the hardest 
one, the successive multiple interactions being then 
ordered by their hardness. There is no initial state 
radiations for interactions in the multiple scattering 
and no color connection between different ladders. 

• For the red closed star symbols, the PYTHIA 6.2 
tune D6 is used. This is the best tune for Tevatron 
data. With the PYTHIA parameter MSTP(82)=4, 
the model for multiple interactions is improved. 
Hadrons are considered to be extended objects. 
The collision has an impact parameter which mod- 
ifies the number of multiple interactions. When the 
protons overlap a lot, there is more chance to have 
multiple interactions than when there is less. 

• The blue closed triangle symbols is PYTHIA 6.4 
all parameters by default, the model for multiple 
interactions is based on the model in PYTHIA 6.2, 
with a modification of the default values for differ- 
ent parameters (hard scale). This is an intermedi- 
ate version between PYTHIA 6.2 D6 and the new 
model. 

• For the purple closed square symbols, a new mul- 
tiple interactions model was developed to be able 
to describe CDF data on {pj) (Fig. [1]). In this 
new model, based on the previous one for the 
determination of the number of interactions, ini- 
tial state radiations are implemented with color re- 
connections between the strings produced by dif- 
ferent interactions. The first interaction is still the 
hardest one. With this new model, more particles 
are produced in event with a lot of multiple interac- 
tions, color re-connections are an important feature 
for the structure of the whole event. 

A detailed description of the physics implied in the 
different PYTHIA models for multiple interactions can 
be found in jstl ■ 

The naive expectation is that, if the production of a 
Jjij/ has nothing to do with multiple interactions, we 
should observe a flat ratio on Fig. |2l It is clearly not 
the case, in all PYTHIA configurations. Furthermore, 
the results depends strongly on the PYTHIA model of 
multiple interactions assumed. Quite unexpected is the 
shape of the strongly decreasing curves we observe, es- 
pecially in the most modern tunes. It is to be noted that 
the multiple interactions scheme in PYTHIA is still un- 
der study and the production of quarkonia is not well 
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Figure 2: PYTHIA Study of Jj^fj production as a function of mul- 
tiplicity, in proton-proton collisions at 10 Tev for different PYTHIA 
settings. For the multiplicity, only charged particles with a pj > 0.9 
GeV/c and |>i < 2.4 are considered. 



constrained, especially with the final state radiation, and 
the production mechanism. 

In any case, the physical meaning (in PYTHIA) of 
the different curves is hard to understand. A detailed 
study of all PYTHIA parameters would be necessary to 
identify what mechanism is responsible in PYTHIA for 
each kind of multiple interactions models, (initial sate 
interaction, final state interaction, pQCD scales, color 
re-connection...). 

More important is to study this observable in real 
data, for instance in the 7 TeV proton-proton data col- 
lected at LHC, and to compare the experimental behav- 
ior to the simulated ones. We do not expect PYTHIA to 
quantitatively predict such a new observable, for which 
the event generator was not built, but rather the experi- 
mental measurement to help tuning PYTHIA. 

To interpret this future measurement, an important 
problem will be its normalization. For corresponding 
study in heavy-ions collisions, one compares to proton- 
proton collisions, scaled by the number of equivalent 
proton-proton collisions. Here we are looking only at 
proton-proton collisions: there is no such underlying el- 
ementary collisions concept. The normalization is then 
not clear and can be debated. 

A possibility would be to compare with open charm. 
Open charm are produced in standard 2 — » 2 hard pro- 
cesses. If J/i// has the same behavior as a standard hard 
processes from the multiplicity point of view, we should 
see a similar behavior. One can also think to look at dif- 
ferent pt classes of the J/ifr to look for eventual modifi- 
cation of the trend looking at hardest particles. Finally, 
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in a more general study, we would recommend to look at 
any observable, one could think of, as a function of mul- 
tiplicity. Taking hadrons or jets for example, one then 
could use another event generator than PYTHIA, with a 
different feature of multiple interactions, and collective 
phenomenon, and here, one could think of EPOS 11311 . 
that does not include, massive quarks and quarkonia 
production. 



4. Conclusion 

This first generator level study of Jlxp production as 
a function of multiplicity shows that the production of 
J/i// in a complete event with the corresponding under- 
lying event and multiple interactions scenario does not 
seem to be under control in the PYTHIA event gen- 
erator Our recommendation is to perform this study 
with the new LHC data at 7 TeV in proton-proton col- 
lisions. The ALICE, ATLAS, CMS and LHCb exper- 
iments have complementary acceptance. For instance, 
CMS can measure J/ip with \y\ < 2.4 but missing the 
lower pt part of the spectrum at midrapidity, while AL- 
ICE will measured J/tfr in rapidity ranges of \y\ < 0.8 
and 2.5 < y < 4, at all pr- This measurement should be 
performed in the four experiments and compared. More 
generally, a detailed study for all observables should be 
performed regarding the multiplicity. 

Recently, the CMS experiment showed an interest- 
ing and unexpected feature of high multiplicity events in 
proton-proton collisions lITsll . The observed long range 
correlations are hard to interpret, but already tells us that 
high multiplicity proton-proton colHsions are not trivial. 
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